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ABSTRACT 


Ultraviolet and optical spectra of in terstellar gas along the li nes of sight to 
nearby stars have been interpreted by iRedheld fc Linskvl (120081) and previons 
stndies as a set of discrete warm, partially ionized clond s each with a d if ferent 
flow vector, temperatnre, and metal depletion. Recently, Grv fc Jenkins fj2014 1 
have proposed a fnndamentally different model consisting of a single clond with 
nonrigid flows hlling space ont to 9 parsecs from the Snn that they propose better 
describes the local ISM. Here we test t hese fnndamenta l ly diff erent morphological 
models against the spatially nnbiased iMalamnt et al.l (1201411 spectroscopic data 
set, and hnd that the mnltiple clond morphology model provides a better £t 
to both the new and old data sets. The detection of three or more velocity 
components along the lines of sight to many nearby stars, the presence of nearby 
scattering screens, the observed thin elongated strnctnres of warm interstellar 
gas, and the likely presence of strong interstellar magnetic helds also snpport the 
mnltiple clond model. The detection and identihcation of interclond gas and the 
measnrement of nentral hydrogen density in clonds beyond the Local Interstellar 
Clond could provide future morphological tests. 
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1. DEVELOPMENT OF MULTIPLE CLOUD MORPHOLOGIES 


The interstellar medium plays a critical role in astrophysics as the interface between 
the creation of metals in stars and their eventual inclusion in the next generation of stars. 
A signihcant fraction of the mass of galaxies is interstell ar gas and the a ccretion and loss 


of this gas play critical roles in galactic evolution (e.g., iKennicuttI Il998l) . Studies of the 


physical processes, thermal structure, and dynamics of interstellar gas include theoretical 
models, numerical simulations, and extensive observational studies typically involving high- 
resolution spectra of ultraviolet absorption lines. The local interstellar medium (LISM) 
which exists in close proximity to the Sun (i.e., within 100 pc), is the closest and simplest 
sample of interstellar gas to study with high-sensitivity absorption line spectra. A signihcant 
advantage of using LISM observations to evaluate interstellar structures is the small number 
of traversed absorbers and high likelihood that the interstellar gas velocity structure along 
the line of sight will be resolved in a high-resolution stellar spectrum. 

From an analysis of h igh-resolu ti on T i II A3384 absorption line spectra toward stars 
within 100 pc of the Sun, ICrutcherl fjl982l) found that warm gas in the LISM is moving 
coherently with a heliocentric velocity vector (/o,6o,K) = (25°,-|-10°, —28 km s“^). He 
noted that this how is consistent with an expanding shell of gas accelerated by the hot stars 
and supernovae in the Sco-Cen Association, and that this interstellar wind how could explain 
the velocity of resonantly scattered solar H I Lyman-a and He I resonance lines. 


Further understanding of the morphology and kinematics of the LISM required more 
sensitive high-resolution spectra, in particular using ultraviolet absorption lines of H I and 
D I Lyman-a, Mg II A2796,2802, Fe II A2600, O I A1302,1304,1306, C II A1334,1335, and 
many weaker transitions of neutral and singly-ionized elements. These spectral lines are 
particularly useful because the ions are abundant in the ISM and the permitted transitions 
are from the highly populated ground states. Such data began to appear in quantity from 
the increasingly sensitive and higher resolution spectrographs on Copernicus, International 
Ultraviolet Explorer (lUE), the Goddard High Resolution Spectrograph (GHRS) instrument 
on the Hubble Space Telescope {HST), and hnally the Space Telescope Imaging Spectrograph 
(STIS) instrument on HST. 


Using Copernicus, iMcGlintock et al.l fjl978[) was able to measure H I column densities 
and radial velocities of interstellar gas in the lines of sight toward eight G and K stars 
within 15 pc of the Sun supporting the model of uniform flow of neutral interstellar gas 
in the s olar neighborh o od. l UE spectra of the H I Lyman-a line toward many hot stars 
allowed iFrisch fc YorkI (119831 1 to create maps of hydrogen column density, which show an 
asymmetrical structure with an V(H I) hole centered at I ~ 225°, b ~ —15°. 
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fiLinskv et al 


_ Analysis of high-resolution GHRS spectra of the nearby st ars a Cen flLinskv fc Wood 

996), Sirius (Lallement et alll994 : Bertin et alJll995 ). Procyon ( Linsky et alll995 ). Capella 


19931) ■ and others proyided accurate column densities, yelocities and nonther- 


mal broadening parameters for these lines of sight by combining information from the sat¬ 
urated H I Lyman-a line, the thermally-broadened but unsaturated D I Lyman-a line, and 
lines of heavy elements like Mg II and Fe II that reveal the velocity structure along the line 
of sight. 

The next step in our understanding of the structure of partially-ionized gas in the LISM 
was the recognition that individual comoving structures (called clouds) could be identihed by 
the common sp ace yelocities of gas acr oss large regions of the sky. Early efforts used Ca II ob¬ 
servations te.g.- ILallement et al.lll986), although these clo u d vectors were signihcantly revised 


with a larger sample of observations fIVallerga et al 


essentially Ca II observations, and iLallement et al. 


I 993 I). Lallement &: Bertin ( 1992 ). using 


(119951) . using primarily high-resolution 


ultraviolet spectra, identified two clouds: the G cloud in the direction of the Galactic Center, 
and the Local Interstellar Cloud (LIC) centered in the opposite di rection and likely c o ntain- 
ing the Sun on the basis of their projected velocity differences. iRedfield fc Linskyl (120001) 
produced a three-dimensional map of the LIC on the basis 16 lines of sight observed by 
GHRS, 13 lines of sight observed with the Ca II lines, and three lines of sight to hot white 
dwarfs observed by the Extreme Ultraviolet Explorer {EUVE). The ac cumulation of add i¬ 
tional observations with th e GHR S and STIS instruments on HS T by iDring et al.l (1l997l ): 
Frisch. Grodnicki. fc Welty ( 2002 ): Redfield fc Linsky ( 2002 . 2004) and o thers were the ba¬ 
sis for more detailed models of the LISM. iFrisch. Grodnicki. fc Weltyi (120021 ) proposed a 
kinematic model with seven clouds based on 96 velocity components detected toward 60 
stars. 


Redfield fc Linskyl (120081 ) [hereafter RL08] next developed a 15 cloud model of the LISM 


on the basis of 270 radial velocity measurements toward 157 stars located within 100 pc of 
the Sun. The assignment of individual velocity components to a cloud was based mainly 
on kinematics (i.e., the radial velocities are consistent to within the measurement errors 
of a common velocity vector for a number of targets distributed over a large region of the 
sky), although contiguity of the sight line coordinates was also a criterion. The shape of the 
clouds, which was drawn by eye to include the G alactic coordinates of the assigned targets, 
is subjective and has been refined over time fe.g.. iMalamut et al.ll2014j [hereafter M14]) with 
the inclusion of new sight lines. Because of the close proximity of these clouds to the Sun, 
they typically subtend large angles on the sky. This provides leverage to measure the three- 
dimensional velocity vector from a sample of radial velocity measurements, and thereby to 
measure the transverse velocities of these clouds. Clouds at a greater distance, or with small 
linear sizes cannot be identified with this technique. The conditions for which a cloud can 
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be identified are complex and hard to qnantify, bnt the accnracy of the velocity vector can 
be estimated from the errors in its components (see Table 16 of RL08). 


All of these warm partially-ionized clonds, now called the complex of local interstellar 
clonds (CLIC), lie entirely or in part within 15 pc of the Snn becanse the nearest stars 
showing absorption by the gas in each clond he within this distance. Since all of these sight 
lines show absorption by partially-ionized gas in these clonds or in not yet identified clonds, 
there is presently no direction in space where logiV(H I) < 17.4 toward an observed star. 



that the inferred temperatures and nonthermal broadening parameters assume a Maxwell- 
Boltzman distribution of velocities and mass-independent turbulence, both of which may not 
be valid in low density clouds. 


Frisch! (120091 1 concluded that the physical properties of the CLIC clouds are typical 


of warm partially-ionized gas observed elsewhere in the solar neighborhood on the basis of 
temperature, velocity, composition, ionization, and magnetic field properties. In particular, 
the ionization equilibr ium of the CLIC gas is consistent with the local EUV radiation field 
fjSlavin fc Frisch! !200^. The recent comprehen s ive re view of the interstellar medium sur¬ 
rounding the Sun by Ihrisch. Redfield. fc Slavin (2011) describes our present understanding 
of the Galactic environment of the LISM, the role of outflowing gas from the Sco-Cen asso¬ 
ciation as a driver for the CLIC kinematics, the interstellar radiation field, ionization and 
depletion of metals in the gas, kinematics of the gas, and the interstellar magnetic field. 


AN ALTERNATIVE MODEL FOR THE CLIC MORPHOLOGY 


Cry fc Jenkins! (120141) [hereafter GJ14] have proposed a fundamentally different mor¬ 


phological model for the warm interstellar gas located near the Sun in which all of the space 
within about 9 parsecs of the Sun consists of a single conti nuous cloud with a nonrig id flow 
and a gradient in metal depletion properties. They fit the !Redfield fc Linskv! (120021 ) Mg II 
and Fe II radial velocity data set with a nonrigid flow that is differentially decelerated in the 
direction of motion and expanding in directions perpendicular to the flow. In their Single Lo¬ 
cal Gloud model, the flow speed and direction near the Sun is consistent with the speed and 
direction of interstellar helium flowing into the heliosphere as measured by the Interstellar 
Boundary Explorer (IBEX) and Ulysses satellites. They found that their single local cloud 
model fits nearly all of the Mg II and Fe II velocity components that RL08 assigned to the 
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Lie, G, NGP, Blue, Leo, Aur, and Get clouds. If the mean neutral hydrogen density is 0.055 
cm“^, as inferred from the hydrogen column densities and distances of nearby stars, then 
their single cloud model hlls all of the space out to roughly 9 parsecs from the Sun. GJ14 
also found that another set of velocity components that RL08 had assigned to the Hyades 
and Mic clouds and including many of the previously unassigned velocity components, could 
be fit by a second vector that they called the Getus Ripple, which may be a signature of a 
shock front inside of the local cloud. They speculated that the remaining unassigned velocity 
components are also perturbations located inside of the local cloud. 


3. A CRITICAL TEST OF THE TWO MORPHOLOGICAL MODELS 


Both the GJ14 and RL08 models were constructed to fit the iRedfield fc Linskvl f 2002h 
data set of Mg II and Fe II absorption line radial velocities, although the RL08 model 
was constructed to also fit radial velocities for other lines of sight measured from Ca II 
line spectra. A critical test of the viability of both models is, therefore, to determine how 
accurately each model predicts the radial velocities of a new data set with targets randomly 
distributed in Galactic coordinates. 


M14 obtained high-resolution STIS E230II spectra of the Mg II, Fe II, and Mn II lines 
in the 2580-2805 A spectral region for stars that previously had only 1200-1700 A spectra. 
This new data set consists of 76 velocity components measured in the lines of sight toward 
34 stars. To be consistent with the RL02 sample, we limit our analysis to stars within 100 
pc, and therefore only used 32 of the 34 sight lines from this sample. The data were obtained 
in the SNAPshot observing mode—designed to provide short (typically one spacecraft orbit) 
observations to fill gaps in the HST schedule. Since the observed targets were selected by the 
HST schedulers from a large list of stars distributed across the sky, the observed targets are, 
in effect, randomly distributed in Galactic coordinates as shown in Figure 1 in M14. The 
M14 data set, therefore, provides us with an appropriate basis for testing the two models. 

We will evaluate how well the GJ14 and the RL08 models £t the new data and utilize 
some statistical tests to make quantitative comparisons, however statistical tests have their 
limitations as the fitting parameters for the two models are very different and could lead to 
possible systematic errors. Figure 1 compares the observed minus predicted radial velocities 
for the two models as a function of Galactic longitude and latitude using the M14 data 
set. The left panels separate by color the predicted velocities for the Local Gloud (32 
measurements) and Getus Ripple (9 measurements) using the procedure proposed by GJ14. 
The right panels make a similar comparison for the exact same set of velocity components, 
except using instead the RL08 model for their suite of discrete clouds. Note that two 
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components associated with the Cetus Ripple by GJ14 are not identihed with a known 
cloud in the RL08 model resulting in 32 + 9 = 41 total points in the GJ14 comparison and 
32 + 7 = 39 points in the RL08 comparison. For the latter comparison, we use a weighted 
mean velocity based on all available LISM measurements, which for the M14 sample is 
primarily Mg II and Fe II. This mitigates systematic errors that can arise from using a single 
ion (e.g., due to saturation or wavelength calibration issues). 


Figure 2 shows a similar set of comp arisons, but now using t he R L08 dataset. This 
datase t includes observations tabulated in iRedheld fc Linskvl fl2002[) and iRedheld fc Linskv 
020041) ■ totaling ~80 sight lines. The GJ14 interpretation of the data results in 80 Local 
Gloud and 36 Getus Ripple assignments, totally 116 components. The RL08 interpretation 
of the data results in 122 cloud identideations. This comparison using this full sample is not 
as clear of a test of the two morphological models because both models were constructed to 
£t these data (together with some Ga II data for the RL08 model). Nevertheless, this second 
comparison is useful as a conhrmation of the results of the hrst comparison. 


As shown in Figures 1 and 2, the RL08 model clearly provides a much tighter £t to 
both the M14 and RL08 data sets. The root mean square (RMS) scatter for the RL08 
model is a factor of 1.5 or two times smaller than for the GJ14 model: 1.4 versus 2.1 for 
the M14 dataset and 0.77 versus 1.5 for the RL08 data set. However, the RL08 model has 
45 free parameters (3 for each of the 15 clouds), versus the 8 free parameters in the GJ14 
model. For the GJ14 model, there are 3 parameters for the Local Gloud velocity vector, 
an additional 3 for the deformation (direction and magnitude), and 2 for the Getus Ripple, 
which is just an offset and range around that offset. An R-test shows that even with more 
free parameters, the RL08 model is statistically preferred. If we look at the entire data set, 
that is data included in RL08 and M14, the GJ14 model has a reduced xt of 49.5, with 149 
degrees of freedom (157 — 8), whereas the same 157 components in the RL08 model has a 
reduced xl of 30.4, with 112 degrees of freedom (157 — 45). There is a 0.4% probability 
that such a dramatic difference in the ratio of could be the result of a random set of data 
as compared to the true model, even taking into account the higher number of parameters 
for the RL08 model flBevington fc Robinsonlll992l) . Note that the x^ calculations above use 
the individual, published velocity errors, whereas GJ14 implemented a velocity error floor of 
1.1 km s“^. If we make the same assumption regarding the velocity errors, the reduced xl 
values become 2.30 and 0.92 for the GJ14 model and RL08 model, respectively. However, 
the basic conclusion remains the same, that the difference in the quality of £t between the 
two models is signiheant. 


Of the 210 velocity components (147 from the RL08 data sample and 63 from the M14 
sample, which excludes those >100 pc and a component identihed with absorption from 
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a circumstellar disk), the RL08 model is able to successfully fit 175 (83%), whereas the 
GJ14 model hts 157 (75%). While the GJ14 model must have a Local Gloud component 
along every sight line, in practice, essentially all sight lines have at least one successfully 
assigned component in the RL08 model. On the other hand, the GJ14 model is unable to 
account for sight lines with more than two absorption components (Local Gloud and Getus 
Ripple), whereas the RL08 model can account for them with several clouds along the line 
of sight. Although the velocity component assignment percentages are similar between the 
two models, the predictive power of the RL08 model is much better than the GJ14 model 
when tested against the M14 data set and the entire data set as a whole. 


For all 34 stars observed and analyzed by M14, the RL08 kinematic model predicts 
40 absorbers along the lines of sight. All 40 radial velocity predictions are successfully 
detected in the M14 sample to within measurement uncertainty (i.e., within ~3cr). The 
RL08 kinematic model also accurately predicts the observed radial velocities for 18 absorbers 
along sight lines lying within 20° of a cloud boundary. The RL08 model clouds have the 
condition of being contiguous, and have a spatial extent (as does the Getus Ripple in the 
GJ14 model). Given the relative sparseness of the sample, the precise cloud morphologies are 
not necessarily well constrained, but will be rehned as new data is made available as long as 
the cloud boundary remains contiguous. In their high-resolution STIS spectra of three early 
B stars located about 70 pc from the Sun, I Welsh fc LallementI fj2010l) also found that 8 out 
of 11 velocity components for these lines of sight are consistent with the predicted velocities 
of six clouds in the RL08 kinematic model. The 40 out of 40 success rate in ascribing M14 
radial velocities to clouds conhrms the robustness of the RL08 model, although any model 
consisting of rigid flow structures is an approximation to what is likely a more complex flow 
pattern that could include rotation in addition to nonrigid translation terms. It is likely that 
a more realistic description of the GLIG kinematics lies somewhere between the two extreme 
morphologies. 


4. DISCUSSION: WHICH MORPHOLOGICAL MODEL IS MORE 
REALISTIC AND WHY IS THIS IMPORTANT? 

Whether the warm partially ionized interstellar gas located within a few parsecs of 
the Sun has a continuous or a discrete multicloud morphology is important because the 
morphology can provide important clues concerning the physical properties of the gas and 
its evolution. 

The presence of discrete clouds requires a structuring agent that could be ionizing radia¬ 
tion from an assorted distribution of different hot stars shielded by clouds in some directions. 
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Fig. 1.— Comparison of the velocity offsets (observed minus model prediction) for the M14 
data set. Left panels use the GJ14 model (red for Component 1 and blue for the Cetus Ripple 
component). Right panels use the RL08 model (the color coding is by cloud structure, which 
is dominated by the LIC [blue], G Cloud [violet], and NCR Cloud [red]). The root mean 
scatter (RMS) is given in the top right corner, and is a factor of 1.5 lower for the RL08 
model than the G14 model. 
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Fig. 2.— Comparison of the velocity offsets (observed minus model prediction) for the 
RL08 data set. Left panels use the GJ14 model (red for Component 1 and blue for the 
Cetus Ripple component). Right panels use the RL08 model (the color coding is by cloud 
structure, which is dominated by the LIC [blue], Hyades Cloud [orange], G Cloud [violet], 
Mic Cloud [green]. Blue Cloud [dark blue], and NGP Cloud [red]). The RMS scatter is given 
in the top right corner, and is a factor of two lower for the RL08 model than the G14 model. 
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inhomogeneous magnetic fields, shock waves, or thermal instabilities. If the warm gas is in¬ 
deed structured into discrete separated clouds, then there must be gas located between the 
clouds that is difficult to detect in spectral lines of hydrogen and singly-ionized metals. This 



An advantage of the continuous cloud morphology model is that it does not require in¬ 
terstitial gas and is therefore consistent with its nondetection, whereas discrete detached 
multicloud models require the presence of such gas and, if/when detected, the properties of 
this gas will provide insight into the energy balance of the CLIC. 


Whether the warm cloud gas fills a small fraction or essentially all of the volume within 
about 9 pc of the Sun depends upon the neutr al hydrogen gas density . If n(III)Ri 0.20 cm“^ 
is typical for the warm gas in nearby clouds, as ISlavin fc Frisch! (1200811 find for the LIC, then 
the warm partially ionized clouds fill a small fraction of nearby space. If, on the other hand, 
n(HI)!^ 0.055 cm“^, as suggested by GJ14, then the warm gas fills all of nearby space and 
there is no interstitial gas within about 9 parsecs of the Sun. 


Gry fc JenkinsI (120141 ) note that if the LISM is comprised of discrete clouds and the 


filling factor is low, it is likely that sight lines could be found that do not traverse any 
warm gas and, therefore, show no LISM absorption. However, every sight line analyzed 


so far shows detected interstellar absorption with logA^(H I) > 17.4. iRedfield fc Linskv 


(120081 ) reported a filling factor of 5.5-19% for their suite of 15 clouds, but this calculation 
assumed the maximum possible spacing of the clouds (limited to the closest stars that showed 
absor ption for a particula r cloud ) and the placement of the clouds within 15 pc of the Sun. 
Also, Welsh fc Lallement ( 201ol) argued that the clouds are all located within 10 pc as 
there is little warm gas between 10 pc and the 70 pc location of their targets. In order to 
explore the question of why there are not gaps in H I absorption, we recomputed the suite of 
simulations described in IRedfield fc Linskvl (120081) while keeping track of the probability of 
any given sight line avoiding all simulated clouds. Our objective was to find an upper limit 
to the total volume enclosing the LISM clouds that predicts that essentially all lines of sight 
to nearby stars have observable HI absorption. We made several simplifying assumptions 
in this analysis and our conclusions are based on a statistical sample of simulated LISM 
configurations. A more accurate estimate would require a full three-dimensional model of 
the LISM based on the entire dataset. We find that if all the clouds are located within 7 pc 






































11 


of the Sun, the hlling factor ranges from 55-100% and there is a ~99.6% probability that 
all sight lines will traverse at least one absorbing cloud (using the original 15 pc limit, this 
probability is ~65%). We conclude that if the 15 clouds in the RL08 model are located within 
7 pc, there should be few or no gaps in the projection of these clouds on the sky and no sight 
lines to nearby stars with undetected H I absorption. Note that this conclusion is strongly 
dependent on the number of very nearby stars (e.g,. <5 pc), and an observational campaign 
to expand this sample could be a powerful discriminant between these two models. While 
both the GJ14 and RL08 models largely hll the immediate interstellar space surrounding 
the Sun, the kinematic properties of the two models are quite distinct, which may hold clues 
regarding the origin and evolution of this material. 


There is evidence that the interstellar ma gnetic held strength immediatel y surround¬ 
ing the heliosphere lies in the range 2.7-5 pG flFrisch. Redheld. fc Slavinll201ll) . Since the 
beginning of 2013, the Voyager 1 spacecraft has been making in situ measurements of the 
interstellar magnetic held strength outside of the termination shock with a mean value of 
4.64±0.09/i G and nearly cons tant orientation very different from that of t he solar magnetic 
held inside of the heliosheath (IGurnett et al.ll2013l: iBurlaea fc Nessll2014aUbl) . Although the 


measured magnetic held likely refers to the held draped around the heliopause, which may be 
compressed relative to the undisturbed held far away from the Sun, Voyager 1 has provided 
us with the best available estimate of the interstellar magnetic held strength near the Sun. 

Since equipartition between magnetic and thermal energy is about 2.7 fiG in the LIG, the 
Voyager 1 magnetic held strength measurement indicates that strong interstellar magnetic 
helds are present near the Sun and could, therefore, control the structure of discrete clouds. 

The thin elongated structures of several clouds identihed by RL08 (Aur, Get, Mic, and 
perhaps others) are consistent with conhnement by elongated parsec scale magnetic helds. 

Large-amplitude intraday and annual scintillations of some well monitor ed, unresolved _ 

uasa rs at radio wavelengths indicate turbulent-scattering screens in the GLIG. lLinskv. Rickett. fc Redheld 
12008 ) showed that the variability of three well-studied quasars (B1257-326, B1519-273, and 
J1819-I-385) can be understood as produced by the Earth’s orbital motion through the dihrac- 
tion pattern of scattering screens located within 7 parsecs of the Sun. The lines of sight to 
these quasars pass close to the outer edges of two or more adjacent or perhaps colliding 
clouds. The shear of the diherent cloud velocities and the likely higher ionization of the 
gas at the cloud edges due to external ionizing radiation or thermal conduction from hot 
surrounding gas could produce the turbulent ionized plasma and the inhomogeneous index 
of refraction properties of scattering screens. The transverse velocities measured for these 
scintillation screens match very well with a cloud in their line of sight for the RL08 model 
but poorly for gas at the Local Standard of Rest. The nonridgid structure of the GJ14 model 
is consistent with a scattering screen in the J1819-I-385 sight line, but cannot explain the 
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scattering screens for the other two sight lines. The existence of nearby scattering screens 
provides additional evidence that isolated warm clouds with ionized edges rather than a sin¬ 
gle cloud even with a nonridgid velocity structure can more naturally explain the complexity 
observed in the LISM velocity structure. 


Another test of whether or not warm partially-ionized gas is confined in identifiable 
clouds, would be the observation of absorption by high ionization species located near the 
outer boundaries of clouds produced by either thermal conduction from surrounding hot gas 
or ionization by the extreme ul traviolet radiation from e CMa and other sources. Searches 
for O VI A1032 absorption (e.g., ISavage fc Lehnerll2006l) led to detections in hot white dwarf 


photospheres but not in the inter stellar medium. I n their summary of searches for lower 
stages of ionization in the LISM, IWelsh et all fj2010jl found that the only clear example is 
interstellar C IV absorption in the line of sight to the B2Ve star HD 158427 located at a 
distance of ~ 74 pc inside of the Local Cavity. This detection is very interesting because 
the observed radial velocity, —24.3 ± 2.0 km s“^, is consistent with the projected radial 
velocities of both the G and Aquila clouds and the line of sight to HD 158427 is tangential 
to the edges of both clouds as shown in RL08. This is the most favorable geometry for the 
detection of weak absorption, because of the long line of sight through the cloud edges. 


5. CONCLUSIONS 

The MI4 data base of interstellar radial velocities provides an unbiased test of the 
robustness of two fundamentally different morphological models of the interstellar medium 
within a few parsecs of the Sun. Both the multiple discrete cloud model proposed by RL08, 
for which each cloud has a rigid flow vector, and the single local cloud model with nonrigid 
flows proposed by GJI4 are approximations to what is likely a more complex kinematic and 
morphological structure. Nevertheless, it is instructive and important to test the predictive 
power of both models. We hnd that the RL08 model hts the new velocity data signihcantly 
better than the GJI4 model and provides a natural way of explaining the observed multiple 
velocity components along the lines of sight to many nearby stars. Also, the likely presence 
of strong magnetic fields and scattering screens in the GLIG are arguments for the presence 
of multiple clouds that do not entirely hll nearby space, although a close packing of clouds 
is likely, resulting in a high filling factor and a low probability of a sight line displaying no 
observable LISM absorption. Important future tests of the GLIG morphology would be the 
detection and identihcation of interstitial gas, if it is present, and the measurement of neutral 
hydrogen densities far away from the heliosphere. A more densely sampled observational data 
set, particularly comprised of very close stars (e.g., <5 pc), would also provide dehnitive 
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tests between smoothly varying deformed kinematics of a continuous medium and the more 
striking variations that would arise in a suite of discrete absorbers. We encourage future 
model development and testing to better understand the physical properties and evolution 
of the interstellar medium close to the solar neighborhood. 
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